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Abstract. Using speckle-interferometry we have carried out repeated measurements of relative positions for the 
components of 34 T Tauri binary systems. The projected separation of these components is low enough that 
orbital motion is expected to be observable within a few years. In most cases orbi tal m otion has indeed been 
detected. The observational data is discussed in a manner similar to Chez et al. (1995). However, we extend 



their study to a larger number of objects and a much longer timespan. The database presented in this paper is 
valuable for future visual orbit determinations. It will yield empirical masses for T Tauri stars that now are only 
poorly known. The available data is however not sufficient to do this at the present time. Instead, we use short 
series of orbital data and statistical distributions of orbital parameters to derive an average system mass that is 
independent of theoretical assumptions about the physics of PMS stars. For our sample this mass is 2.0 Mq and 
thus in the order of magnitude one expects for the mass sum of two T Tauri stars. It is also comparable to mass 
estimates obtained for the same systems using theoretical PMS evolutionary models. 
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1. Introduction 

The mass is the most fundamental parameter of a star, 
because to a large extent it determines its structure and 
evolution. Therefore, it is a major problem for our under- 
standing of pre-main sequence (PMS) evolution that at 
this time there are no reliable empirical mass determina- 
tions for young stars with M < 1 Mq. Mass estimates for 
T Tauri stars are usually based on comparisons of their 
position in the Hertzsprung-Russell diagram (HRD) with 
theoretical PMS evolutionary models, which means that 
they are affected by the (unknown) uncertainties within 
these models. Moreover, it is not possible to rate the qual- 
ity of different sets of PMS models by comparison to ob- 
servational data. 

For this reason, empirical mass determinations for young 
stars are highly desirable. Binary stars offer a unique pos- 
sibility to do this, because the system mass is known as 
soon as the orbital parameters are determined. There are 
many binaries among T Tauri stars in nearby star-forming 
regions (SFRs). Most of them have been detected during 



the last decade by high-angular resolution surveys in the 
near infrared (NIR) (For an overview of this topic, see the 
review by Mathieu et al. ( 200C| ) and references therein.). 
The first reliable empirical masses of PMS stars were given 
by Casey et al. ( 1998) ) for the components of the eclipsing 
double-lined spectroscopic binary (ESB2) TYCrA. These 
masses are M 1 = 3.16±0.02 M Q and M 2 = 1.64±0.01 M . 
The secondary mass is consistent with th e pred ictions of 
PMS models fr om D 'Antona & Mazzitelli ( |l994| ) and also 
Swenson et al. ( 1994 ). The primary is already close to the 
main sequence. The lowest-mass PMS stars with empir- 
ically determined masses thus far known are the com- 
ponents of RXJ0529.4+0041. For this ESB2, Covino et 
al. Q200CD determined the masses Mi = 1.25±0.05 Mq and 
Mi = 0.91 ± 0.05 Mq. They concluded that these masses 



to: 



Jens Woitas, e-mail 



Send offprint requests 
woitas@tls-tautenburg . de 

* Based on observations collected at the German-Spanish 
Astronomical Center on Calar Alto, Spain, and at the 
European Southern Observatory, La Silla, Chile. 



are in good agreement with with the Baraffe et al. (1998) 
and Swenson et al. (1994) models, but less consistent with 
sets of PMS tracks from D'Antona & Mazzitelli fll994| ) 
and Palla & Stahler (1999). Because of the relatively high 
masses, these results cannot be used to check the PMS 
models for K- or M-dwarfs and objects with masses below 
the hydrogen burning mass limit at w 0.075 Mq. 
In this paper we will follow the approach of Ghez et 
al. ( 1995| , herafter G95). Using NIR speckle interferome- 
try, they obtained repeated measurements for the relative 
astrometry of the components in 20 T Tauri binary sys- 
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Fig. 1. The first row shows modulus (left) and 
bispectrum-phase (right) of the complex visibility for the 
binary XZ Tau, derived from data obtained on 29 Sep 1996 
at the 3.5 m-telescope on Calar Alto with the NIR array 
camera MAGIC at A = 2.2 fim. The second row is the 
modulus and phase of an artificial complex visibility that 
is fitted to the data. This fit is indicated in the third row in 
a one-dimensional projection towards the connection line 
of the components (perpendicular to the strip patterns in 
Fourier space) . The circle around the strip patterns corre- 
sponds to the Nyquist spatial frequency, 7.0arcsec _1 for 
the adopted pixel scale. 



terns. In this way they showed that in most of these sys- 
tems, orbital motion can be determined. From short pieces 
of orbital data and the statistical distribution of orbital 
parameters, they have derived the average system mass of 
1.7 Mq, which is in the order of magnitude expected for 
the mass of two T Tauri stars. 

We present similar data for 34 T Tauri binary systems and 
in this way increase the object list and also the observa- 
tional time base. An overview of our observations and data 
reduction is given in Sect. H. The results are presented in 
Sect.f| discussed in Sect.^fand summarized in Sect.||. 



2. Observations and data reduction 



2.1. The sample 

We have repeatedly observed 21 system s in T aurus- Auriga 
detected as binaries by Leinert et al. (1993) and Kohler 
& Leinert (1998) and 11 systems in Scorpius-Centaurus 
detected by Kohler et al. (2000). Furthermore, we re- 
observed two binaries found by Ghez et al. (1997a). These 
are HM Anon in the ChamaeleonI association, and HN 
Lup in the Lupus SFR. We combine our data with results 
taken from literature (see Sect. |2.6| ). In particular twelve 
systems discussed by G95 are also objects of our study. 



2.2. Observations 

The objects in Taurus- Auriga were observed with the 
3.5 m-telescope on Calar Alto. After September 1993 
these observations used the NIR array camera MAGIC. 
Previous measurements were obtained with a device 
for one-dimensional speckle interferometry described by 
Leinert & Haas ( 1989Q . The observations of young bina- 
ries in southern SFRs were carried out at the ESO New 
Technology Telescope (NTT) at La Silla that is also a 
3.5 m-telescope. The instrument used for these observa- 
tions was the NIR array camera SHARP I of the Max- 
Planck Institute for Extraterrestrial Physics (Hofmann et 
al. 1992 ). Both cameras are capable of obtaining fast se- 



quences of short exposures with integration times r « 
0.1s, which is cru cial for the applied data reduction pro- 
cess (see Sect. 2J3). 

Most of the data were obtained in the K-band at 
A = 2.2 fim. Some observations used the J-band at 
A = 1.25 /zm and the H-band at A = 1.65 /im. In 
these cases the main goal of the observations was to obtain 
resolved photometry of the components at those wave- 
lcnghts. In the course of data reduction we could however 
show that highly precise relative astrometry can also be 
derived from observations in J and H (see Sect. 2.4 for the 
determination of binary parameters). 

2.3. Speckle interferometry 

At the presumed dis tance of the objects (between 140 and 
190 pc, see Sect. 3T and Tablejl]) orbital motion measur- 
able within a few years can only be expected for the clos- 
est pairs with projected separations of d < 0'.'5. Therefore, 
high angular resolution techniques are necessary that over- 
come the effects of atmospheric turbulence and yield the 
diffraction-limited information about the objects. For this 
purpose, we have used speckle interferometry. 
Sequences of w 1000 short exposures (r sa 0.1s) are 
taken for the object and a nearby point source, the ref- 
erence star. The integration time is shorter than the co- 
herence time of the turbulent atmosphere, so the turbu- 
lence is "frozen" , and the images are noisy, but principally 
diffraction limited. After Fourier transforming these "data 
cubes", the power spectrum of the image is deconvolved 
with that of the reference star to obtain the modulus of 
the complex visibility. The phase is reconstructed using 
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the Knox-Thompson algorithm (Knox & Thompson 1974) 
and the bispectrum method (Lohmann et al. 1983). The 
complex visibility is the Fourier transform of the object 
brightness distribution. For a sufficiently bright object it 
will contain the diffraction-limited information. 



2.4. Determination of binary parameters 

Modulus and phase of the complex visibility are charac- 
teristic strip patterns for a binary, as shown in Fig. |] (first 
row) for the XZ Tau system. The position angle of the 
companion is orientated perpendicular to these patterns 
and towards higher values of phase. The distance between 
two strips is inversely proportional to the projected sepa- 
ration of the components, and the amplitude of the pat- 
terns denotes the flux ratio. Position angle, projected sep- 
aration and flux ratio are determined by constructing an 
artificial complex visibility from a set of these parame- 
ters (second row in Fig.[l]) and fitting it to the data (third 
row in Fig.0). This fit uses the Amoeba algorithm (Press 
et al. 1994) . The errors of the binary parameters are esti- 



Table 1. Adopted distances to stars in nearby SFRs 



mated by applying this fitting procedure to different sub- 
sets of the data. 



2.5. Pixel scale and detector orientation 

To obtain position angle and projected separation as as- 
tronomical quantities they must be transformed from the 
array onto the sky by calibrating them with pixel scale 
and orientation of the detector. Minimizing the uncer- 
tainties that result from this calibration is crucial for a 
highly accurate determination of the components' rela- 
tive astrometry that is proposed here. We determine pixel 
scale and detector orientation from astrometric fits to im- 
ages of the Orion Trapezium cluster where precise astrom- 
etry has been derived by McCaughrean & Stauffer ( 1994 ). 
Typical errors are 0.1° for the detector orientation and 
10~ 4 arcsec/pixel for the pixel scale that is « O'.'07/pixel 
for MAGIC and w O'.'05/pixel for SHARP I. 
This precise calibration exists for all observations obtained 
by the authors since July 1995. For calibrating data from 
previous observing runs we used visual binary stars with 
well-known orbits. We observed these stars afterwards and 
calibrated their separation and position angle with the 
help of the Trapezium cluster. Thus, we have placed all 
our (two-dimensional) speckle observations into a consis- 
tent system of pixel scale and detector orientation. 



2.6. The database 



In Table |A.1| the derived position angles, projected sepa- 
rations and flux ratios are presented. For the subsequent 
analysis (Sect.||) the relative positions are combined with 
data taken from literature. Twelve of our systems have 
also been discussed by G95. For some binaries there are 
additional measurements obtained by HST imaging, the 
HST Fine Guidance Sensors or adaptive optics (see refer- 



SFR 



Distance [pc] Reference 



Taurus-Auriga 142 ± 14 

Scorpius-Centaurus 145 ± 2 

Chamaeleonl 160 ± 17 

Lupus 190 ± 27 



1998 



Wichmann et al (|1998 
de Zeeuw et al. (|199S|) 
Wichmann et al 
Wichmann et al 



1998 



ences m Table [Al]). The finding that our relative astrom- 
etry fits very well with that obtained by other groups fur- 
ther supports the idea that the position angles and pro- 
jected separations are as precise as is indicated by the 



errors given in Table Al 



3. Results 

3.1. Relative motion of the components 



In Fig. A.l the relative positions of the components at 
different epochs are shown in cartesian and polar coordi- 
nates. These plots are only given for the 23 out of our 34 
systems for which there are at least three data points. To 
obtain a simple approximation of the relative velocity we 
applied weighted linear fits to this data. For the 11 sys- 
tems with only two observations we simply connect the 
two data points and derive the error of the relative ve- 
locity from the uncertainties of the two separations and 
position angles. 

For a quantitative analysis, the relative velocities must 
be transformed to an absolute length scale. This requires 
knowing the distances of the discussed objects. We adopt 
distances to the SFRs that are the mean of all Hipparcos 
distances derived for members of the respective associa- 
tion. The values and references are given in Tabled There 
remains, however, an uncertainty, because distances of in- 
dividual objects may be different from these mean values. 
To take this into account, we assume that the radial di- 
ameters of the SFRs are as large as their projected di- 
ameters on the sky. The latter quantity can be estimated 
to be ~ 20° for Taurus- Auriga (see Fig. 1 in Kohler & 
Leinert |1998) as w ell as Scorpius-Centaurus (see Fig. 1 in 
Kohler et al. 2000 ) . Concerning the mean distances from 
Table 0, this corresponds to a diameter of 50 pc. Thus we 
will assume ±25 pc as the uncertainty for the distance of 
an individual system, which is an upper limit: More than 
two thirds of the stars will be within ± 15 pc for an even 
distribution. 

The velocities of the companions relative to the primaries 
derived by applying the assumed distances are given in 
Table H They are also plotted in Fig.|| in cartesian and 
polar coordinates (similar to Fig. 3 in G95). Our measure- 
ments can only cover the projection of motion onto the 
sky, so the v p are given with respect to the main com- 
ponent, not to the observer. The adopted v is the mean 
of the respective values derived from the fits in cartesian 
and polar coordinates. In 3 out of 34 systems v is different 
from zero on the 3<r level, in 9 systems on the 2tr level and 
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Table 3. Distribution of the companions with respect to 
certain relative velocities and separations 





d < 25 AU 


d > 25 AU 


v > 5km/s 


7 ±2.6 


11 ±3.3 
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3 ±1.7 


13 ±3.6 
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Fig. 3. Relative velocities of the components as function 
of their mean projected separation. 

in 18 systems on the la level. Thus, we are fairly confident 
that there really is relative motion of the components in 
most systems. 

3.2. Origin of the relative motion 

We now examine the origin of this relative motion. For 
this purpose, we must discriminate orbital motion from an 
apparent relative motion that can be caused by the proper 
motion of a T Tauri star with respect to a background 
star projected by chance or by the proper motions of two 
T Tauri stars projected by chance. One has further to 
consider that "companions" to T Tauri stars detected with 
only one observation in one filter are not necessarily stellar 
and may be Herbig-Haro objects. We will also examine the 
possible influence of unresolved additional components on 
the observed motion. 

3.2.1. Orbital motion 

First we will derive an upper limit for relative velocities 
caused by orbital motion. That limit is given by the con- 
dition that in the case of orbital motion the kinetic energy 
of the system is less than its (negative) potential energy 
and is equal to it in the extreme case of a parabolic orbit, 
i. e. 



2 2GMs ys 
v < — . 



(1) 



We adopt Mtts < 2 Mq as upper mass limit for one 
T Tauri star (Hartmann 1998) and take the mean of the 
measured projected separations as estimate for r. There 
is only one companion with a relative velocity that is 
larger than the value derived from Eq.^j, namely that of 
RX J1546. 1-2804. The relative velocity of this companion 
is, however, still consistent with orbital motion, consider- 
ing the la error in its v. The lower limit for v is zero be- 
cause in short pieces of orbit as discussed here the orbital 
motion may occur purely radial to the observer. Thus, the 
large majority of the velocities from Tableware consistent 
with orbital motion. 

Furthermore, it is interesting to examine the relationship 
between separation and relative velocity. In the special 
case of a circular orbit observed face-on this relation will 
be 



V (X 



(2) 



where r is the instantaneous separation of the compo- 
nents. 



for the orbital velocity v and the semimajor axis a. 
Eq.|| will be modified by projection effects and the actual 
orbital parameters, and the unknown system mass is an 
additional parameter in this relation. However, close com- 
panions should, on average, move faster than distant ones. 
In Fig. U the measured projected velocities, v, are plotted 
as a function of the components' projected separation, d. 
The correlation is weak for the reasons mentioned above, 
but there is at least a tendency to fulfil the prediction of 
Eq.|^. Among the closest companions with d < 25 AU, rel- 
ative velocities with v > 5km/s are more frequent than 
for the companions with larger separations (Table|^). 
We conclude from this section that the observed relative 
velocities in almost all cases are not in contradiction to 
orbital motion. For a final classification, other possible 
origins of the relative motion must be considered. 

3.2.2. Chance projected background stars 

A background star projected by chance will probably 
be located at a much larger distance than the observed 
T Tauri star, so its proper motion can be neglected. The 
measured relative motion is thus expected to be the proper 
motion of the T Tauri star. In Taurus- Auriga there are 
subgroups with different directions of proper moti on an d 
a mean proper motion of 15.8 km/s (Jones & Herbig 1979). 
For Scorpius-Centaurus, de Zeeuw et al. ( 1999 ) give a 
mean proper motion of 17.6 km/s with a mean direction of 
v x = -17.4 km/s and v y = -2.8 km/s. For both SFRs the 
distances given in Table were adopted. There are only 
two systems with relative velocities in this order of mag- 
nitude, namely V 773 Tau (v = 13.78 ± 7.34 km/s) and 
RX J1546.1-2804 (v = 23.60 ± 6.28 km/s). 
For V 773 Tau, the Hipparcos catalogue gives proper mo- 
tions of fi x — —24.89 ± 1.89mas/yr and fi y = 0.65 ± 
2.83mas/yr. The resultant v — 16.68 ± 2.28 km/s is com- 
parable with the observed relative velocity of the com- 
ponents, but the proper motion of V 773 Tau happens 
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Table 2. Projected relative velocities of the companions with respect to the primaries in cartesian and polar coordi- 
nates. The adopted v is the mean of the total velocities derived from (v x , v y ) and (v p , v$). Note that v p is given with 
respect to the main component and not relative to the observer. 



System 


«i[km/s] 


v y \km/s\ 


v p [km/s] 


»4W S 1 


v[km/s 


] 


d[AU] 


V 773 Tau 


-1.35 ± 2.47 


15.64 ± 6.87 


-8.37 ± 6.82 


8.41 ± 2.83 


13.78 ± 7.34 


14 


9 


LkCa 3 


4.54 ± 1.92 


1.16 ± 1.07 


-0.60 ± 1.49 


-2.98 ± 1.62 


3.86 


± 


2.20 


68 


7 


FO Tau 


3.67 ± 1.63 


-5.62 ± 1.52 


-1.48 ± 1.36 


5.63 ± 1.72 


6.27 


± 


2.21 


22 


3 


CZ Tau 


-4.83 ± 1.46 


-1.97 ± 2.22 


-2.72 ± 2.27 


4.98 ± 1.43 


5.45 


± 


2.67 


45 


3 


FS Tau 


-9.26 ± 2.87 


2.06 ± 1.44 


-1.88 ± 1.54 


6.61 ± 2.77 


8.18 


± 


3.19 


35 


6 


FW Tau 


8.11 ± 2.12 


-7.69 ± 2.22 


-7.77 ± 2.24 


6.87 ± 1.67 


10.77 ± 2.93 


15 


8 


LkHa 331 


-3.77 ± 2.21 


2.51 ± 1.04 


-1.55 ± 1.65 


1.76 ± 1.77 


3.44 


± 


2.43 


40 


3 


XZ Tau 


2.14 ± 1.03 


4.25 ± 1.32 


-0.35 ± 0.50 


-5.03 ± 1.52 


4.90 


± 


1.64 


43 


2 


HK Tau G2 


0.41 ± 0.34 


-2.06 ± 1.92 


5.00 ± 1.80 


-2.05 ± 1.92 


3.75 


± 


1.95 


26 


3 


GG Tau Aa 


-0.62 ± 0.69 


-6.59 ± 1.82 


0.41 ± 0.72 


-4.71 ± 1.80 


5.67 


± 


1.94 


35 


8 


UZ Tau/w 


3.20 ± 1.31 


1.36 ± 1.24 


1.21 ± 1.32 


2.76 ± 1.24 


3.24 


± 


1.81 


50 


9 


GH Tau 


8.61 ± 2.96 


-2.97 ± 2.14 


-1.42 ± 3.18 


-3.19 ± 1.92 


6.30 


± 


3.68 


45 





Elias 12 


-13.79 ± 3.95 


-0.65 ± 1.71 


-6.05 ± 3.67 


-4.68 ± 2.56 


10.73 ± 4.39 


49 


3 


IS Tau 


-4.55 ± 2.35 


0.45 ± 1.92 


0.21 ± 1.93 


3.24 ± 2.37 


3.91 


± 


3.04 


31 


7 


IW Tau 


-2.50 ± 0.55 


-1.14 ± 1.91 


2.50 ± 0.47 


2.79 ± 1.80 


3.24 


± 


1.97 


39 


6 


LkHa 332/G2 


2.31 ± 1.80 


7.68 ± 3.23 


-5.58 ± 2.46 


3.72 ± 2.07 


7.36 


± 


3.45 


36 


5 


LkHa 332/G1 


-1.66 ± 1.62 


0.54 ± 1.23 


2.86 ± 1.06 


5.50 ± 1.76 


3.97 


± 


2.04 


32 


2 


LkHa 332 


2.91 ± 3.54 


1.21 ± 1.20 


-0.30 ± 3.69 


-0.83 ± 0.60 


2.02 


± 


3.74 


47 





BD+26 718B Aa 


-6.98 ± 2.14 


5.40 ± 1.66 


-8.26 ± 2.70 


-0.35 ± 0.12 


8.55 


± 


2.71 


67 


5 


BD+26 718B Bb 


0.98 ± 3.60 


1.06 ± 1.51 


2.07 ± 3.48 


3.02 ± 1.84 


2.55 


± 


3.92 


23 


3 


BD+17 724B 


0.49 ± 1.64 


10.83 ± 6.20 


-0.41 ± 3.88 


-4.02 ± 5.15 


7.44 


± 


6.43 


12 


8 


NTTS 155808-2219 


0.75 ± 2.64 


-4.06 ± 2.06 


3.43 ± 2.13 


2.31 ± 1.72 


4.13 


± 


3.05 


29 


4 


NTTS 155913-2233 


1.60 ± 0.80 


-5.38 ± 2.03 


3.09 ± 0.90 


-5.30 ± 1.85 


5.87 


± 


2.11 


43 


3 


NTTS 160735-1857 


0.89 ± 2.51 


-6.38 ± 2.24 


0.17 ± 1.91 


-6.44 ± 2.28 


6.44 


± 


3.18 


43 


4 


NTTS 160946-1851 


-1.02 ± 0.64 


0.09 ± 0.61 


0.18 ± 0.65 


1.04 ± 0.64 


1.04 


± 


0.90 


30 


5 


HM Anon 


1.04 ± 9.33 


-7.04 ± 8.94 


-2.65 ± 7.43 


6.61 ± 6.42 


7.12 


± 


11.37 


42 


1 


HN Lup 


7.42 ± 3.30 


-0.94 ± 3.21 


1.57 ± 3.33 


-7.32 ± 3.15 


7.48 


± 


4.60 


46 


3 


RXJ1546. 1-2804 


22.34 ± 1.58 


2.86 ± 4.86 


-3.77 ± 1.84 


-24.40 ± 7.21 


23.60 ± 6.28 


14 


3 


RXJ1549.3-2600 


-0.83 ± 1.33 


-0.45 ± 1.23 


-0.17 ± 1.40 


0.92 ± 0.46 


0.94 


± 


1.64 


23 


7 


RXJ1600.5-2027 


-4.02 ± 1.61 


1.10 ± 1.42 


1.88 ± 1.69 


-3.72 ± 1.21 


4.17 


± 


2.12 


28 


2 


RXJ1601. 7-2049 


2.04 ± 1.75 


1.56 ± 1.73 


0.00 ± 1.37 


-2.57 ± 1.77 


2.57 


± 


2.11 


29 


7 


RXJ1601.8-2445 


8.43 ± 4.55 


1.46 ± 4.76 


7.79 ± 3.53 


-3.59 ± 5.74 


8.56 


± 


6.66 


13 


3 


RXJ1603.9-2031B 


-2.03 ± 2.93 


-9.58 ± 1.73 


-4.28 ± 2.11 


-8.95 ± 3.19 


9.86 


± 


3.61 


15 


8 


RXJ1604.3-2130B 


-2.98 ± 1.82 


-2.52 ± 1.43 


1.03 ± 1.55 


3.78 ± 1.67 


3.91 


± 


2.13 


12 


3 



almost only in declination (X in Fig. A.l ) which is contra- 
dictory to our observations. Furthermore, for V 773 Tau 
the relative motion is still explainable with orbital motion, 



considering the upper limit derived in Sect. 3.2.1 



In the case of RXJ1546. 1-2804 there are no Hipparcos data 
for this individual object, so we adopt the mean proper 
motion for the OB association Upper Scorpius given by 
de Zeeuw et al. ( 1999| ). These values are \i = 17.7km/s, 
[i x — —17.4 km/s and fi y — —2.8 km/s with (formal) er- 
rors of pa 0.1 km/s. The respective values from Tablc|] 
are v = 23.6 ± 6.3 km/s, v x = 22.3 ± 1.6 km/s and 
v y = 2.9 ±4.9 km/s. This is a close correspondance, given 
that the the proper motion presented by de Zeeuw et 



al. (1999) is not for this single object, but for the whole 
association. Note that fi and v must be antiparallel if the 
companion is a chance-projected background star, because 
in that case fi is the motion of the primary with respect 
to the "companion" . Another argument that supports the 
idea that the companion of RX J1546. 1-2804 is a chance- 
projected background star is that its relative motion is 



above the limit given by Eq.[l]. It is, however, still consis- 
tent with orbital motion at the la level. 
One must take into account that in both cases the mea- 
sured projected separations are d ~ O'.'l, whi ch ma kes any 
chance projections very unlikely (see Sect. 3.2.3 for the 
probability of chance projections). Thus, one must con- 
sider other origins of these high relative velocities. One 
solution may be that the distances of individual systems 
are largely different from the adopted values (Table[l]) for 
the SFRs as a whole. The unusually high derived values for 
relative velocities would then be caused by overestimat- 
ing the objects' distances. Another possibile explanation, 
namely the presence of unresolved additional companions, 



will be discussed in Sect. 3.2.5 



The companion of RX J1546. 1-2804 may be such a chance- 
projected background star. In the case of V 773 Tau this 
seems to be unlikely because the direction of proper mo- 
tion does not match, however, the high relative velocity 
of the components remains problematic. Further obser- 
vations of these systems will be necessary to determine 
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Fig. 2. Relative velocities of the components in T Tauri binary systems in cartesian coordinates (left panel) and polar 
coordinates (right panel). The dashed line in the right panel points towards the locus of Herbig-Haro objects that is 
far outside this plot. RX J1546. 1-2804 is not plotted in this figure. Its locus is out of the picture in both panels. This 
plot is similar to Fig. 3 in G95. 



whether there is a curvature in the relative motion which 
would undoubtedly classify it as orbital motion. In gen- 
eral, chance-projected background stars are not frequent 
among close visual companions to T Tauri stars. 

3.2.3. Chance projection of two T Tauri stars 

There may also be companions that are in fact objects 
projected by chance belonging to the same SFR as the 
"primary" . In that case both components would have 
roughly the same proper motion. Any relative motion 
would then be caused by the velocity dispersion within the 
respective SFR or subgroup, as was mentioned by G95. 
These velocity dispersions are between 1 and 2 km/s in 
each coordinate for d ifferent subgroups in Taurus- Auriga 
(Jones fc H crbig 197S ). For Scorpius-Centaurus, de Zeeuw 
et al. ( 1999 ) formally derived a velocity dispersion of only 
O.lkm/s. One cannot distinguish this kind of relative mo- 
tion from slow orbital motion by analyzing relative veloc- 
ities. 

It is, however, not probable that chance projections within 
the same SFR are a frequent phenomenon, because of the 
low stellar density in the SFRs discussed here. Leinert et 
al. (1993) concluded that there are less than 4-10 -5 objects 
brighter than K = 12 mag per arcsec 2 in Taurus- Auriga. 
This includes association members and field stars. All 
companions discussed here are brighter than K = 12 mag 
and have projected separations of less than 0'.'5. The mean 
number of chance-projected companions within a radius of 
0'.'5 around our 21 objects in Taurus- Auriga is thus 



jV proj = 21 • 0.5 2 • 7T • 4 • 1(T 5 = 6.6 • 1CT 4 . 



(3) 



For Scorpius-Centaurus, Kohler et al. ( |200C ) derived a 
stellar density of (6.64 ± 0.45) • 10~ 4 arcsec -2 from star 
counts in the vicinity of X-ray selected T Tauri stars. 
Applying Eq.|| yields a mean number of 5.7 • 10 -3 chance- 
projected companions around 11 objects in Scorpius- 
Centaurus. These values are low enough to rule out that 
chance projections within the same SFR cause the detec- 
tion of an appreciable number of close "companions" . The 
low observed relative velocities in some systems are prob- 
ably due to the fact that we can only detect the projected 
motion on the sky. If most of the orbital motion happens 
radial to the observer in the time span covered by our data 
it will not be measurable by our observations. 



3.2.4. Herbig-Haro objects 

G95 raised the question of whether an appreciable num- 
ber of companions of young stars detected with only one 
observation in one broad band filter may in fact not be 
stellar, but rather condensations in gaseous nebulae. Such 
companions would appear as Herbig-Haro objects that are 
driven by strong winds of active T Tauri stars. They usu- 
ally have velocities of some 100 km/s (e. g. Schwartz 1983). 
Moreover, their relative motion with respect to their driv- 
ing source can be only radial, not tangential. Thus, the lo- 
cus of Herbig-Haro objects in Fig. || (right panel) is in the 
direction of the dashed line in this figure, but far outside 
the plot. As was reported by G95, we find no system where 
the observed relative velocity is consistent with these re- 
strictions. 

Furthermore, in the case of HV Tau, where companion 
C has been declared as a Herbig-Haro object, we could 
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show that in fact it is a stellar companion (Woitas & 
Leinert 



1998) 



3.2.5. Apparent relative motion caused by unresolved 
companions 

Another problem with observed relative velocities in mul- 
tiple systems has been pointed out by G95. If the main 
component has an additional unresolved companion, or- 
bital motion in this close pair can shift the photocen- 
ter of the "primary" . This will be misinterpreted as mo- 
tion of the visual secondary if only relative astrometry 
is measured, as is the case here. G95 have used this ef- 
fect to explain the surprisingly high relative velocity in 
the Elias 12 system where the visual companion is mov- 
ing with 10.73 km/s at a separation of 49.3 AU (Table|J). 
This is only slightly below the upper limit for orbital mo- 
tion, which is 12.0 km/s in this case (Eq.|l|). In this system, 
the primary has another companion with a separation of 
d = 23 mas, detected by Simon et al. ( |1995 ) using lu- 
nar occultations. If the system mass of the close pair is 



ctcrs. The results presented in Table [A~l| , however, remain 
valuable for future orbit determinations that will yield em- 
pirical masses for T Tauri binary systems. Furthermore, 
it is already possible to estimate an average system mass 
from this database. This average mass is not dependent 
on theoretical assumptions about the physics of PMS evo- 
lution and should therefore be a reliable empirical estima- 
tion of the masses of T Tauri stars. To derive this mass 
we follow the approach of G95, but improve it in some 
important aspects. 

First, we write Kepler's third law in the natural units 
M/Mq, a/1 AU and P/lyr: 



MSy 



a" 

P2 



(4) 



and assume a face-on circular orbit where the total orbital 
motion happens tangential to the observer and therefore 
equals the observed velocity: 



^face— on.circ 



2na 

~P~' 



assumed to be 1M and a relation < d >= 0.95 a be- In this special case (i = 0°, e = 0) Eq.| becomes 



tween the mean projected se parat ion < d > and the semi- 
major axis a (Leinert et al. 1993| ) is adopted, one would 
expect the period of the close pair to be P « 6.4 yr. This 
is roughly the sa me t imespan that is covered by our ob- 
servations (Table AT). The shift of the photocenter of the 
close pair with respect to the visual secondary will cancel, 
and thus the derived secondary's velocity is not affected 
by the tertiary. 

Also in the V 773 Tau system, where we have noticed an 
unusually high relative velocity of the visual secondary 
(Sect. fT2~2) the re is an additional spectroscopic compan- 
ion (Wclty 1995 ). The period of this close pair is 51.1 days, 
so possible shifts of the photocenter are less than 2 mas 
and thus not measurable by our observations. 
A candidate for a system where the observed relative 
motion may be influenced by an unresolved tertiary is 
BD+26 718B Aa. In this system v — 8.55 km/s at a sep- 
aration of 67.5 AU. Similar to Elias 12, this value is close 
to the upper limit for orbital motion (v < 10.2 km/s from 
Eq. 0) , but also far below the relative velocity expected for 
a background star projected by chance (Sect. 3.2.2). 



Msy S = 



face- 



4tt 2 



(5) 



(6) 



Unfortunately, we don't know a or i>f a cc-on.circ of our 
systems. G95 used computer simulations of binary orbits 
to derive statistical relations between the observed quan- 
tities dp r oi and v h s on the one hand and a and Vf t 



acc — on,circ 



on the other hand. They assumed an eccentricity distri- 
bution of /(e) = 2e and an isotropic distribution of incli- 
nations and obtained the following results: 



< d pro j >= 0.91a . 

< «obs >= 0.72u face _ c 



(7) 
(8) 



3.2.6. Conclusions 

The observed relative velocties can be explained by or- 
bital motion. Together with the result that other origins 
of relative motion can be ruled out at a high confidence 
level for most systems, we conclude that the observed 
motion is orbital motion in nearly all binaries discussed. 



3.3. Estimation of an empirical average system mass 

For all binary systems discussed here the available por- 
tions of the orbit are too short to calculate orbital param- 



However, d pro j and i> bs are not uncorrelated. For ex- 
ample, if we observe a companion in the outer part of its 
eccentric orbit, d pro j will be larger, and v bs smaller than 
their average values. Therefore, we cannot insert these 
equations into Eq.|| to obtain a statistical relation be- 
tween dp ro j, «obs, and M Sys . 

To overcome this problem, we performed more sophis- 
ticated computer simulations. Each simulation contains 
10 million binaries with a fixed system mass and ran- 
domly distributed orbital parameters. The periods fol- 
low the distribution of periods of main-sequence stars 



(Duquennoy & Mayor 1991), the distribution of eccen- 
tricities is /(e) = 2e and the inclinations are distributed 
isotropically, while all the other parameters have uniform 
distributions. The distances to the observer are varied 
within a range of 143±25pc. We chose two observation 
dates separated by a random timespan between 4 and 
10 years and computed the average projected separation 
and orbital velocity, in much the same way as we did for 
the real data. We then select binaries in the projected 
separation range from 10 to 70 AU. For these binaries, we 
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compute < v 2 -d > /Ms ys . These simulations are repeated 
for different system masses. For Ms ys between 0.5 Mq and 
2.5 M , the results vary from 18.4 to 18 AU 3 /(yr 2 M ). 
This gives us following relation: 



M 



1 



Sys 



18.2 



< v 2 ■ d > 



(9) 



The term on the right hand side of Eq.^ contains v 2 . 
Since the measured velocities and separations differ from 
the real values by some unknown measurement error, this 
leads to an additional bias term: 



< v lbs 



> = < (u rca i + Sv) 2 ■ (d rca i + Sd) > 

= < V 2 cal ■ d rca i > + < 5v 2 ■ d rea l > ■ (10) 



Combining Eq.|| and 10, and using Av as an estimate 
for Sv, we obtain a relation that allows us to use the ob- 
served values to obtain an average system mass: 

1 ' ' ' (11) 



Ms 



18.2 



\<v, 



obs ' dproj 



> 



< A «obs • dproj 



> 



This is only an estimate for the mass in an ensem- 
ble of systems with identical system masses. The stars 
in our sample will probably not all have the same mass. 
However, we expect them to have similar masses, and for a 
sufficiently large number of systems, the uncertainties will 
cancel. The result should yield a reliable average mass for 
these systems. 

If we use the data of all stars in our sample, Eq. [TT| 
yields a system mass of 2.5 Mq. Excluding RXJ1546.1- 
2804 lowers the result to 2.0 Mq. If we further exclude 
V 773 Tau, Elias 12, and BD+26718BAa, we arrive at a 
system mass of 1.3 Mq. We have reason to assume that 
the companion to RXJ1546. 1-2804 is a chance projected 
background star. For the three other systems mentioned, 
the velocity is puzzling at first sight, but is still consistent 
with orbital motion. Furthermore, other possi ble ex pla- 
natio ns fail for V 773 Tau and Elias 12(Sects. ^2~^ and 
3.2.5 ). Thus, it does not seem justified to exclude them 
from the sample, and we will adopt 2 Mq as the result for 
the average dynamical system mass. 

Given the statistical uncertainties, it is difficult to esti- 
mate the error of the system mass. Using the standard 
deviation of the quantities averaged in Eq 
the error of the mean yields ss 0.7 Mq. This is in agree- 
ment with the scatter we obtain if we exclude or include 
the stars mentioned in the last paragraph. 
We conclude that the average mass of the systems in our 
sample is in the range 1.3. . . 2.5 Mq, with a most probable 
system mass of 2 Mq . Our result is thus consistent with 
the expectation that T Tauri stars' masses are around 
M w 1 Mq and also with the average mass of 1.7 Mq 
that G95 derived for their sample. 



4. Discussion 

4.1. Comparision with theoretical results 

We have estimated masses for the components by compar- 
ison with theoretical PMS models for a subsample that 



contains 17 out of our 34 systems. In these cases, we ob- 
tained resolved J band photometry. This spectral band 
is supposed to be least affected by circumstellar excess 
emission and can thus be taken as an indicator of stel- 
lar luminosity. Moreover, in these 17 systems, there are 
no additional companions known, so the mass sum of the 
components derived from theoretical models will match 
the dynamical system masses. We placed the components 
of these systems into the HRD, estimating the stellar lumi- 
nosity from the resolved J band magnitudes, assigning the 
optical sp ectral type of the system (taken from Keny on fc 
Hartmann 1995j for Taurus- Auriga and Walter et al. 1994 
for Scorpius-Centaurus) to the primary and assuming that 
all components within one system are coeval. Masses of the 
components were then derived using t he PM S evolution- 
ary tracks of D'Antona & Mazzitclli ( 1995 ), Swenson et 
al. (|1994|) and Baraffe et al. fll998| ). 

The mean mass obtained for this subsample fro m the 
HRD is 0.88 M for the DAntona & Ma zzitell i fll998| ) 
tracks and 1.17 Mq for the Swenson et al . ( 1994 ) tracks. 
The PMS model from Baraffe et al. ( 1998 ) yields a mean 
mass of 1.28 Mq. The uncertainties of the mass estimates 
that originate from observational data used for placing 
the components into the HRD are w 0.1 Mq. The aver- 
age empirical mass derived from Eq.[ll] for that subsam- 
pleQ is 1.22 ± 0.50 Mq. Within the large formal error, the 
predictions of all three PMS models match our empirical 
result. However, our dynamical < M > is much closer to 



the mean masses derived using the Baraffe et al. (1998) 
and Swenson et al. (1994) tracks than to the value calcu- 
lated from the D'Antona & Mazzitelli fll998j ) PMS model, 
which seems to underestimate T Tauri star masses. This 
fin ding h as been recognized by other authors. Bonnell et 
al. ( 1998[) estimated T Tauri stars' masses from infall ve- 
locities of accreted material and also conclude that the 
empirical masses are generally l arger than those predicted 
by the D'Antona & Mazzitelli (|l998|) mo del. A similar re- 
sult was reported by Simon et al. ( 2000 ) , who calculated 
T Tauri stars' masses from Keplerian motion in circum- 
stellar and circumbinary disks. 



to estimate 4.2. Implications for binary statistics 



There is significant relative motion in most systems, and 
this motion is in almost all cases consistent with orbital 
motion (Sect.||). As already pointed out by G95, this 
demonstrates that the large majority of all close com- 
panions detected in the multiplicity surveys mentioned 
in Sect. 2A really are gravitationally-bound stars. No bi- 
nary component discussed here h as to be reclassified as 
a Herbig-Haro object (Sect. 3.2.4 ), and only 2 out of 34 
companions may be chance projected background stars 
(Sect. |3.2.2l ). 



For the reasons mentioned above, the subsample does 
not contain RX J1546.1-2804, V 773 Tau, Elias 12 and 
BD+26718B Aa. Thus, the mean m ass for the subsample is 
lower than the value derived in Sect. 3.3. 
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This finding is particularly important because there is a 
companion overabundance among T Tauri stars in the 
SFRs discussed here when compared to main-sequence 
stars in the solar neighbourhood (Leinert et al 
Ghez et al. 11993], Ghez et al.|l997al Kohler & Leinert 



1993 



1998 



Kohler et al. pOOC ). In Taurus- Auriga, almost all T Tauri 



stars seem to be components of multiple systems. To fur- 



ther confirm this result, Kohler & Leinert (1998) per- 
formed stellar counts in the vicinity of their survey ob- 
jects and concluded that in Taurus-Auriga, statistically 
4.3 out of 44 apparent comp anion s are projected back- 
ground objects. Kohler et al. (200C), in a similar way, de- 
rived a number of 7.8 chance projections per 46 com- 
panions in Scorpius-Centaurus. We found one candidate 
for a chance-projected background star out of 21 objects 
in Taurus- Auriga and one candidate out of 11 objects in 
Scorpius-Centaurus. The results are not directly compa- 
rable because we have only studied the closest pairs for 
which chance projections are least probable. The percent- 
age of background stars projected by chance among the 
observed companions is, however, of the same order of 
magnitude in both studies. Thus, it can be concluded that 
chance projections do not affect the binary statistics sig- 
nificantly in the SFRs discussed here. 
Kohler et al. (2000) excluded six close companions in 
Scorpius-Centaurus from a restricted sample. Their ob- 
served separations are less than the strict diffraction limit 
X/D = 0'/13 of a 3.5 m-telescope in the K-band, so they 
cannot be definitely distinguished from elongated single 
objects. In three of these cases, namely RX J1601. 8-2445, 
RX J1603.9-2031B and RX J1604.3-2130B we derive a rel- 
ative velocity that is consistent with orbital motion (see 
Sect.||). Thus, we propose to classify these objects as bi- 
nary systems in further studies of multiplicity in the OB 
association, Scorpius-Centaurus. 

5. Summary 

Based on repeated measurements of the relative astrom- 
etry in 34 close T Tauri binary systems, we have repro- 
duced the results given by G95 and extended their work 
to a larger number of binaries and particularly to a much 
longer timespan of up to 10 years. We showed that in most 
systems significant relative motion of the components has 
occured. In almost all cases this relative motion can be 
explained by orbital motion. In only two systems the ob- 
served motion may be the result of the proper motion of 
a T Tauri star that is accidentally projected in the close 
vicinity of a background star. 

From the short pieces of orbit available at the moment 
(up to 20° in position angle), we derive a mean dynam- 
ical system mass of 2.0 ± 0.7 Mq for our sample. This 
mass is consistent with the predictions of current sets of 
PMS evolutionary models within the uncertainties. The 
large formal error of this mean mass does not allow a sig- 
nificant discrimination between different models, but we 
draw the tentative conclusion that t he masses predicted 
by the D'Antona & Mazzitelli ( 1998 ) model may be sys- 



tematically too low. 

The result that orbital motion can be detected in most 
systems discussed here indicates that the "companions" 
found in previous multiplicity surveys really are gravi- 
tationally bound stars. This is a further confirmation of 
the binary overabundance in Taurus- Auriga and Scorpius- 
Centaurus compared to nearby main sequence stars. 
Furthermore, the detection of orbital motion allows a def- 
inite classification of three objects with very close separa- 
tions as stellar companions. 
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Appendix A: Overview of observations 



Table A.l presents an overview of all measurements of the 
relative astromctry in T Tauri binary systems that form 
the database for the present study. In Fig. AT these mea- 
surements are plotted as a function of epoch for the indi- 
vidual systems in cartesian and polar coordinates. 
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Table A.l. Relative astrometry of the components in young binary systems measured at different epochs 
reference is given the results were obtained by the authors. 



System 



Date 


Filter 


Flux ratio J2//1 


PA/° 


dproj / mas 


Reference 


3 Oct 1990 


K 




295 ± 4 


112 ± 1 


G95 


21 Sep 1991 


K 


0.13 ± 0.04 


295 ± 3 


170 ± 10 




10 Oct 1992 


K 




307 ± 5 


120 ± 20 


G95 


5 Oct 1993 


J 


0.113 ± 0.005 


308.9 ± 1.4 


111 ± 4 




25 Nov 1993 


K 




304 ± 8 


92 ± 8 


G95 


19 Oct 1994 


K 




318 ± 2 


65 ± 1 


G95 


29 Oct 1994 


UBVRI (HST) 




321.9 ± 2.7 


63 ± 2 


Ghez et al 


6 Dec 1990 


K 


0.5 


78 ± 1 


470 


± 40 




5 Oct 1993 


J 


0.62 ± 0.02 


74.5 ± 0.7 


494 


± 13 




29 Sep 1996 


H 


0.885 ± 0.011 


74.2 ± 0.3 


485 


± 4 




19 Nov 1997 


K 


0.755 ± 0.007 


73.4 ± 0.1 


485 


± 4 




19 Sep 1991 


K 


0.92 ± 0.04 


180 ± 4 


165 


± 5 




18 Oct 1991 


K 




181.7 ± 0.6 


161 


± 1 


G95 


5 Oct 1993 


J 


0.30 ± 0.01 


182.0 ± 0.9 


183 


± 5 




20 Oct 1994 


K 




190.6 ± 0.4 


153 


± 2 


G95 


14 Dec 1994 


K 


0.72 ± 0.02 


189.7 ± 0.9 


159 


± 4 




19 Dec 1994 


K 




191.2 ± 0.4 


154 


± 2 


G95 


9 Oct 1995 


K 


0.648 ± 0.018 


193.1 ± 1.3 


156 


± 4 




27 Sep 1996 


K 


0.628 ± 0.019 


194.3 ± 0.9 


152 


± 4 




29 Nov 1996 


J 


0.55 ± 0.05 


200.0 ±1.1 


143 


± 4 




16 Nov 1997 


H 


0.698 ± 0.014 


198.9 ± 0.8 


150 


± 5 




19 Mar 1992 


K 


0.46 ± 0.03 


84 ± 3 


330 


± 10 




28 Sep 1996 


K 


0.183 ± 0.004 


90.3 ± 0.3 


306 


± 4 




29 Nov 1996 


J 


0.120 ± 0.005 


89.7 ± 1.1 


320 


±6 




20 Sep 1989 


K 




60 ± 5 


265 


± 5 


Simon et a 


25 Jan 1996 


RI (HST) 




84.0 ± 1.5 


239 


± 5 


Krist et al 


28 Sep 1996 


H 


0.183 ± 0.008 


84.4 ± 1.6 


238 


± 4 




29 Nov 1996 


J 


0.188 ± 0.007 


83.3 ± 3.0 


265 


± 15 




17 Nov 1997 


K 


0.138 ± 0.005 


83.3 ± 1.5 


248 


± 5 




17 Oct 1989 


K 


1.00 ± 0.01 


158 ± 2 


151± 1 




31 Mar 1990 


K 




160 ± 5 


160± 20 


Simon et a 


13 Dec 1994 


K 


0.61 ± 0.10 


187.4 ± 1.6 


98 ± 6 




9 Oct 1995 


K 


1.00 ± 0.05 


188.5 ± 3.7 


78 ± 5 




27 Sep 1996 


H 


0.76 ± 0.10 


190.1 ± 3.3 


71 ± 4 




29 Oct 1991 


K 


0.73 ± 0.04 


290 ± 4 


300 


± 10 




26 Jan 1994 


J 


0.706 ± 0.023 


286.2 ± 0.9 


279 


± 4 




9 Oct 1995 


K 


0.66 ± 0.03 


286.6 ± 0.8 


279 


± 4 




29 Sep 1996 


H 


0.70 ± 0.1 


287.8 ± 0.2 


278 


± 4 




17 Oct 1989 


K 


0.35 ± 0.03 


154 ± 3 


300 


± 20 




10 Oct 1991 


K 




151 ± 2 


310 


± 10 


G95 


27 Jan 1994 


J 


1.51 ± 0.03 


146.1 ± 0.6 


307 


± 4 




28 Jan 1994 


K 


0.41 ± 0.01 


146.8 ± 1.0 


306 


± 5 




19 Dec 1994 


K 




147.0 ± 0.4 


296 


± 2 


G95 


5 Jan 1995 


RI (HST) 




147.8 ± 0.25 


307 


± 9 


Krist et al 


29 Sep 1996 


K 




145.4 ± 0.3 


299 


± 4 




30 Nov 1996 


J 


3.54 ± 0.31 


145.7 ± 0.7 


309 


± 10 




19 Nov 1997 


K 


0.316 ± 0.007 


140.9 ± 0.5 


302 


± 4 




27 Sep 1991 


K 


0.88 ± 0.03 


300 ± 4 


190 


± 10 




19 Oct 1991 


K 




304.1 ± 0.5 


163 


± 5 


G95 


26 Jan 1994 


J 


0.764 ± 0.023 


302.5 ± 0.5 


176 


± 4 




27 Jan 1994 


H 


0.85 ± 0.12 


301.5 ± 1.7 


180 


± 10 




8 Dec 1994 


V 




302.4 ± 0.6 


178 


± 2 


Simon et a 


19 Dec 1994 


K 




301.2 ± 0.6 


181 


± 2 


G95 


9 Oct 1995 


K 


0.85 ± 0.06 


300.3 ± 1.0 


188 


± 4 




27 Sep 1996 


H 


0.76 ± 0.05 


299.6 ± 0.6 


193 


± 4 




19 Nov 1997 


K 


0.587 ± 0.017 


298.3 ± 0.6 


214 


± 4 





V 773 Tau 



LkCa3 



FO Tau 



CZTau 



FS Tau 



FWTau 



LkHa 331 



XZTau 



HK Tau G2 
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System 



Date 


Filter 


Flux ratio I2/I1 


PA/° 


dproj / mas 


2 Nov 1990 


K 


0.64 ± 0.1 


9 ± 2 


255 


± 10 


21 Oct 1991 


K 




2 ± 1 


260 


± 10 


23 Dec 1994 


IJHK 




3 ± 2 


260 


± 10 


27 Jan 1994 


J 


0.543 ± 0.004 


357.8 ± 0.4 


246 


± 4 


25 Jul 1994 


UBVRI (HST) 




358.8 ± 0.5 


250 


± 2 


24 Sep 1994 


K 




357 ± 2 


258 


± 4 


18 Oct 1994 


K 




0.9 ± 0.5 


242 


± 3 


8 Oct 1995 


K 


0.54 ± 0.02 


356.9 ± 0.7 


247 


± 4 


29 Sep 1996 


H 


0.556 ± 0.007 


355.5 ± 0.4 


245 


± 4 


16 Nov 1997 


K 


0.564 ± 0.004 


353.6 ± 0.1 


247 


± 5 


10 Oct 1998 


K 


0.476 ± 0.005 


350.7 ± 0.4 


260 


± 4 


8 Oct 1990 


K 




± 8 


340 


± 60 


10 Nov 1990 


K 




357.4 ± 0.3 


359 


± 2 


26 Jan 1994 


J 


0.76 ± 0.07 


1.3 ± 1.0 


358 


± 10 


24 Jul 1994 


UBVRI (HST) 




0.1 ± 0.1 


368 


± 1 


19 Dec 1994 


K 




1.3 ± 0.4 


360 


± 3 


29 Sep 1996 


K 


0.66 ± 0.02 


1.3 ± 0.2 


366 


± 4 


19 Oct 1991 


K 




118.4 ± 0.1 


314 


± 1 


27 Oct 1991 


K 


0.91 ± 0.04 


120 ± 1 


350 


± 10 


27 Jan 1994 


J 


0.89 ± 0.01 


116.4 ± 0.3 


313 


± 3 


19 Dec 1994 


K 




116.0 ± 0.4 


307 


± 2 


29 Sep 1996 


H 


1.03 ± 0.10 


111.4 ± 1.8 


300 


± 10 


10 Nov 1990 


K 




332.5 ± 0.2 


371 


± 1 


27 Oct 1991 


K 


0.45 ± 0.02 


333.9 ± 0.8 


411 


± 10 


5 Jan 1993 


H 


0.54 ± 0.01 


332.4 ± 0.8 


340 


± 10 


26 Nov 1993 


K 




334 ± 6 


340 


± 20 


27 Jan 1994 


J 


0.579 ± 0.022 


328.5 ± 0.9 


347 


± 5 


27 Aug 1994 


V 




328.5 ± 0.4 


335 


5 ± 2 


24 Sep 1994 


K 




326.9 ± 0.9 


336 


± 3 


18 Oct 1994 


K 




328.1 ± 0.6 


329 


± 4 


4 Feb 1995 


V 




328.8 ± 0.4 


337 


5 ± 2 


29 Sep 1996 


H 


0.58 ± 0.04 


325.2 ± 0.4 


318 


± 4 


19 Oct 1991 


K 




92 ± 2 


221 


± 4 


10 Jan 1993 


K 


0.21 ± 0.02 


86 ± 8 


210 


± 20 


26 Jan 1994 


J 


0.32 ± 0.01 


93.5 ± 0.5 


238 


± 4 


9 Oct 1995 


K 


0.165 ± 0.004 


96.3 ± 1.0 


235 


± 4 


28 Sep 1996 


H 


0.20 ± 0.01 


97.0 ± 0.2 


223 


± 4 


29 Nov 1996 


J 


0.21 ± 0.01 


95.3 ± 1.5 


210 


± 9 


28 Oct 1991 


K 


0.91 ± 0.04 


177 ± 2 


270 


± 20 


27 Jan 1994 


J 


1.24 ± 0.04 


176.0 ± 1.0 


278 


± 5 


28 Sep 1996 


H 


1.30 ± 0.15 


179.4 ± 0.7 


288 


± 7 


27 Oct 1991 


K 


0.60 ± 0.05 


243 ± 2 


300 


± 10 


27 Jan 1994 


J 


0.458 ± 0.003 


236.4 ± 0.3 


260 


± 4 


15 Dec 1994 


K 


0.530 ± 0.017 


237.5 ± 1.0 


252 


± 4 


28 Sep 1996 


H 


0.411 ± 0.012 


240.1 ± 0.7 


239 


± 4 


17 Nov 1997 


K 


0.465 ± 0.022 


241.6 ± 0.4 


234 


± 4 


4 Oct 1990 


K 




77.6 ± 0.4 


208 


± 2 


27 Oct 1991 


K 


0.58 ± 0.03 


85 ± 2 


230 


± 20 


25 Nov 1993 


K 




83 ± 2 


230 


± 10 


18 Oct 1994 


K 




82.7 ± 0.4 


224 


± 2 


12 Dec 1994 


K 


0.555 ± 0.006 


83.8 ± 0.3 


228 


± 4 


28 Sep 1996 


H 


0.560 ± 0.009 


86.3 ± 0.4 


228 


± 4 


29 Nov 1996 


J 


0.511 ± 0.026 


87.4 ± 0.3 


236 


± 5 


22 Nov 1997 


K 


0.51 ± 0.02 


91.0 ± 0.3 


230 


± 4 


27 Oct 1991 


K 


0.23 ± 0.03 


204 ± 2 


360 


± 10 


27 Jan 1994 


J 


0.80 ± 0.04 


204.5 ± 1.0 


313 


± 11 


28 Sep 1996 


H 


0.485 ± 0.015 


204.8 ± 1.3 


325 


± 10 


22 Nov 1997 


K 


0.46 ± 0.01 


203.5 ± 0.3 


343 


± 4 



Reference 



GG Tau 



UZTau/v 



GHTau 



Elias 12 



IS Tau 



IWTau 



LkHa332/G2 



LkHa332/Gl 



LkHa332 



G95 

Roddier et al. (|l996|) 



Ghez et al. ( |l997b| ) 

G95 
G95 



Simon et al. fll992| ) 

G95 



Ghez et al. ( |l997b| ) 

G95 

G95 



G95 
G95 

G95 



Simon et al. fll996| ) 

G95 
G95 



Simon et al. (1996) 



Ghez et al. (1993) 



G95 

G95 
G95 
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System 



Date 



Filter Flux ratio h/h PA/° 



dproj / mas Reference 



BD+26 718B - Aa 



BD+26 718B - B 



BD+17 724 



HM Anon 



HNLup 



16 Sep 1994 
27 Sep 1996 
18 Nov 1997 
16 Sep 1994 
27 Sep 1996 
18 Nov 1997 
8 Oct 1995 
27 Sep 1996 
18 Nov 1997 

4 May 1994 

5 May 1998 

27 Apr 1994 
7 May 1998 



K 
K 
K 
K 
K 
K 
K 
K 
H 
K 
H 
.1 
K 
H 
J 



320.1 ± 0.5 



0.178 ± 0.004 



0.310 ± 0.006 



0.240 ± 0.006 
0.205 ± 0.004 

0.628 ± 0.013 
0.831 ± 0.044 



320.0 
319.9 
136.8 
133.4 
135.2 
208.3 
193. 
190.7 
236 ± 5 
243.6 ± 



0.1 
0.4 
9.0 
0.2 
0.2 
9.7 
2.1 
1.8 

1.0 



359 ± 2 
351.3 ± 0.2 



496 ± 3 
474 ± 3 
454 ± 5 
166 ± 7 
155 ± 3 
171 ± 4 
98 ± 11 
84 ± 3 
87 ± 5 
270 ± 30 
256 ± 7 

240 ± 10 
247 ± 4 



Kohler & Leinert 
Kohler & Leinert 

Kohler & Leinert 
Kohler & Leinert 

Kohler & Leinert 
Kohler & Leinert 



Ghez et al. (1997a 



Ghez et al. (1997a 



RXJ1546. 1-2804 


1 May 1994 


K 




41.2 ± 0.5 


110 ± 3 


Kohler et al 




7 May 1998 


H 


0.537 ± 0.016 


318.7 ± 0.5 


88 ± 4 




RXJ1549.3-2600 


1 May 1994 


K 




320.3 ± 0.4 


164 ± 4 


Kohler et al 




7 May 1998 


H 


0.443 ± 0.003 


322.2 ± 0.2 


163 ± 4 




RXJ1600.5-2027 


2 May 1994 


K 




171.7 ± 0.5 


189 ± 4 


Kohler et al 




7 May 1998 


H 


0.735 ± 0.021 


165.3 ± 0.3 


200 ± 4 




RXJ1601. 7-2049 


2 May 1994 


K 




324.7 ± 0.9 


205 ± 4 


Kohler et al 




7 May 1998 


H 


0.504 ± 0.012 


320.5 ± 0.2 


205 ± 4 




RXJ1601. 8-2445 


10 Jul 1995 


K 




289.6 ± 10 


76 ± 5 


Kohler et al 




7 May 1998 


H 


0.595 ± 0.030 


280.4 ± 1.5 


108 ± 4 




RXJ1603.9-2031B 


2 May 1994 


K 




140.9 ± 0.6 


121 ± 4 


Kohler et al 




7 May 1998 


H 


0.666 ± 0.022 


113.3 ± 1.3 


96 ± 4 




RXJ1604.3-2130B 


1 May 1994 


K 




327.4 ± 0.8 


82 ± 4 


Kohler et al 




7 May 1998 


H 


0.471 ± 0.004 


342.3 ± 0.7 


88 ± 4 




NTTS 155808-2219 


2 May 1994 


K 




313.7 ± 1.2 


193 ± 5 


Kohler et al 




7 May 1998 


H 
J 


0.575 ± 0.008 
0.585 ± 0.019 


317.5 ± 0.8 


213 ± 4 




NTTS 160735-1857 


2 May 1994 


K 




84.1 ± 0.3 


299 ± 4 


Kohler et al 




8 May 1994 


H 
J 


0.584 ± 0.021 
0.656 ± 0.019 


76.9 ± 1.0 


300 ± 7 




NTTS 155913-2233 


8 Jul 1990 


K 




347.7 ±0.7 


283 ± 4 


G95 




4 Jul 1993 


K 




345 ± 2 


292 ± 8 


G95 




1 May 1994 


K 




345.6 ± 0.4 


301 ± 4 


Kohler et al 




7 May 1998 


H 
J 


0.577 ± 0.016 
0.548 ± 0.014 


337.6 ± 0.4 


318 ± 4 




NTTS 160946-1851 


5 Aug 1990 


K 




160.4 ± 0.6 


208 ± 2 


G95 




4 May 1991 


K 




162 ± 3 


212 ± 5 


G95 




12 May 1992 


K 




160 ± 4 


210 ± 20 


G95 




25 Jul 1992 


K 




158.6 ± 0.5 


212 ± 5 


G95 




1 Jul 1993 


K 




159 ± 1 


209 ± 2 


G95 




2 May 1994 


K 




161.9 ± 0.4 


203 ± 6 


Kohler et al 




8 May 1998 


H 


0.194 ± 0.009 


161.4 ± 1.5 


220 ± 13 





al. ( |2000| ; 



al. ( |2000| ; 



al. fl2000[ 



al. ( |2000| ; 



al. ( |200Ct 



al. ( [2000] ; 



1998) 



1998 ) 



1998) 



1998) 



1998) 



1998) 



0.210 ± 0.007 
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Fig. A.l. Relative astrometry of the components in T Tauri binary systems in cartesian coordinates (first and second 
column) and polar coordinates (third and fourth column). The solid lines indicate the results of weighted linear fits to 
this data. Triangles denote our new data presented in this paper. Asterisks refer to measurements of "first epoch" that 
coincide with the detection of the companions. They have already been published (Leinert et al. 1992 and Kohler & 



Leinert 1998). Squares indicate data points taken from literature, in most cases from G95 (see Table A.l for reference) 
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